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ABSTRACT 

We use the Mitchell Spectrograph (formerly VIRUS-P) on the McDonald Observatory 2.7m Harlan J. Smith 
Telescope to search for the chemical signatures of massive elliptical galaxy assembly. The Mitchell Spectrograph 
is an integral-field spectrograph with a uniquely wide field of view (107 x 107 sq arcsec), allowing us to achieve 
remarkably high signal-to-noise ratios of ^20-70 per pixel in radial bins of 2-2.5 times the effective radii of the 
eight galaxies in our sample. Focusing on a sample of massive elliptical galaxies with stellar velocity dispersions 
(T*> 150 km s"', we study the radial dependence in the equivalent widths (EW) of key metal absorption lines. By 
twice the effective radius, the Mgb EWs have dropped by ^ 50%, and only a weak correlation between ct, and 
Mgb EW remains. The Mgb EWs at large radii are comparable to those seen in the centers of elliptical galaxies 
that are ^ an order of magnitude less massive. We find that the well-known metallicity gradients often observed 
within an effective radius continue smoothly to 2.5Re, while the abundance ratio gradients remain flat. Much like 
the halo of the Milky Way, the stellar halos of our galaxies have low metallicities and high a-abundance ratios, 
as expected for very old stars formed in small stellar systems. Our observations support a picture in which the 
outer parts of massive elliptical galaxies are built by the accretion of much smaller systems whose star formation 
history was truncated at early times. 



1. INTRODUCTION 

Elliptical galaxies are comprised of mostly old stars, contain 
little gas or dust, and show very tight scaling relations between 
their sizes, central surface brightnesses, and stellar velocity 
dispersions (the Fundamental Plane, e.g., iDiorgovski & DavisI 
11987; Dressier et al. 1987). Despite their apparent simplicity, 
observations of elliptical galaxies continue to surprise and con- 
found us. Their central stellar populations suggest that the 
most massi ve elliptical galaxies formed their stars rapidly and 
'Fabeij 1 19731 



early (e.g., 



van Dokkum et all2008 . 
2008t iDamianov et all l2009i: iCassata et al.1 120ia) . 



iThomas etani2005l) . And yet, evi- 
dence for dramatic size evolution has emerged, such that ellip- 
tical galaxies at z w 1 were apparently a f actor of ^ 2 smaller 
at fixed mass than they_are today (e.g., 'Truiillo et a l.l l2006t 
\^an der Wei et al. 2008; Cima tti et all 

It is, of 

course, extre mely challenging to measure galaxy size s at high 
redshift (e.g.. lHopkins et al.ll2009HSaracco et al.ll2010l) . but ev- 
idence continues to mount that the size evolution is real (e.g ., 
iNewman et al ] |201 iHBrodie et al.ll201 UlPapovich et al.ll201 ll) . 

The most common scenario to explain the dramatic size 
growth in elliptical galaxies at late times invokes minor merg- 
ing that ca n make galaxies fluffier without adding very much 
mass (e.g., Gallagher & Ostrikei^'1972'; 'Bovlan-K ojchin & Mai 
12007; Naabetal.'2007, 2009; Newman etal. 201 Naively, 
late-time merging with small systems would wash out the well- 
established scaling relations between stellar velocity dispersion 
(cr,) and stellar population properties observed i n local ellip- 
tical g alaxies, such as the Mg/j-cr, relation (e.g., ' Bender et al.l 
Il993h Furthermore, stellar population studies of local ellipti- 
cal galaxies clearly find that the stars in the most massive ellip- 
tical galaxies were formed earliest (z > 2) and fastest (< Gyr), 
while lower-mass systems have mor e extended formatio n histo- 
ries and later formation times (e.g.. iThomas et al. I I2005I) . From 



the tight color-magnitude relation alone it is hard to support 
much late-time star formation (or the addition of more metal- 
poor stars, e.g.. lBower et al ] [T99l . 

The tension between the tight scaling relations of elliptical 
galaxies and their apparent puffing up from late-time merg- 
ing is alleviated if the stars added at late times are deposited 
at large radius. The vast majority of stellar population work 
is heavily weighted towards the very luminous central compo- 
nent of these gal axies, usually well within the half-lig ht radius 
{Re\ e.g.. Faberl 



iKuntschner et al 



19731: iRawle et al.1 l2008t iGraves et al. 2009; 



201( 



Recently, tha nks to the Sloan Digi- 
^n^Aet al."2000^, 

w available for exan unmg (. 

Zibetti et al.' 2005; T ortora etal] 120101; ISuh et al.1 



tal Sky Survey (SDSS;[%Aet al."2000^, very large samples of 
elliptical gal axies are now a vailable for exaimnmg color gra- 
dients (e.£ 



2010t iG onzalez -Perez et al. 201 1) . With few exceptions (e.g., 
Rudick et al. 201^ these observations have not extended much 



beyond the effective radius. To fully exploit the fossil record to 
understand the assembly of elliptical galaxies, we ought to look 
for radial changes in the stellar population, particularly beyond 

Re- 

The study of the radial dependence of chemical composition 
in elliptical galaxies has a long h istory. I maging studies of ellip- 
tical galaxy colors date back to Ide Vauco uleurs (1961). Since 
then, there have been many studie s made of the radial color 
gradients in elliptical galaxies (e.g.. lTiff3ll969l; IWirth & Shawl 
1983; Eisenhardt et al. 2007). The summary presented in 
IStrom & S troml(ll97 8) remains accurate today; elliptical galax- 
ies are bluer at large radii, most likely du e to a decline in metal- 
Ucity jSpinradll 19721; IStrom et al.|[T976l) . However, with pho- 
tometry alone it is difficult to precisely disentangle the well- 
know n degeneracies between age and metallicity (e.g., Wor thevI 
|1994|) . Many spectroscopic surveys have looked at the gradients 
in the equivalent widths ( EW) of key metal lines that can break 
these degeneracies (e.g., ISpinrad & TavloJ [19711: iFaber et aTl 
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19771; iGorgas et alJ ' 1 9901: iFisher et al J 1 1 9951: iMehlert et al 




Broughet alJ l2007t" 



Baes et al 
— 



2003 : [Ogandoetal 

200 7*: Annibalietal. 2007; Sanchez-Blazquez et al.' '20a 
Rawle et al. 2008; Kuntschner et al. 2010; Weiimans et al 
20091) . Roughly speaking, spectroscopic work confirms the 



sured by the SDSS pipeline) and redshifts in the narrow range 
Table 1 . The Sample 



overall conclusions from imaging studies. Metallicity domi- 
nates the color changes, decreasing outwards by 0.1-0.5 dex 
per decade in radius. In general there is no strong evidence for 
age gradients (although see Baes et al. for an alternate view). 

Spectroscopic surveys can track more than just metallicity 
and age. They can also study the relative abundances of in- 
dividual elements. In particular, the a elements (e.g., Mg, C, 
O, N) are formed in Type 11 supernova explosions, while the 
Fe-peak elements (Fe, Cr, Mn) are formed predominantly in 
Type la supernovae, and are thus produced with a temporal lag 
from the peak of star formation. The relative quantity of a 
to Fe-peak elements provides a star-formation timescale, with 
enhanced a/Fe ratios pointing to rapid time-scales of star for- 
mation. Elliptical galaxies display a strong trend of increasing 
a/Fe abunda nce with increasing mass or stellar velocity dis- 
persio n (e.g., lFaberf ll973': 'Terlevich etal."1981|; IWorthev etal] 
Il99l although see also Kelson et al. (200Q. It is therefore 
thought that the most massive elhptical galaxies fo rmed their 
stars rapidly and at z > 2 (e.g.. iThomas et af]|2005b . Thus far, 
no strong g radients in a/Fe ratios have been detected at large 
radii (e.g., Kuntschner et alJUoiOl: ISpolaor et al and ref- 

erences therein). 

A few studies have managed to probe stellar populations in 
elliptical galaxies beyond the effective radius. It is very hard to 
achieve the required signal-to-noise at large radii with long-slit 
spectroscopy since the area subtended on the sky is small and 
the sky level is fac t ors of several brighter than the signal (e.g., 
iKelson et all 120021: ISanchez-Blazquez et all l2007h . Integral- 
field unit (IFU) spectroscopy provides two-dimensional infor- 
mation, and coadding the signal in annuli strongly boosts the 
signal relative to the sky. A handful of studies thus far have used 
IFUs with smaller fields of view, and either tile the instrument 
at large radius (e.g., Weijmans et al. 20 091) or focus on the cen- 
tral re gions of the galaxy ( Rawle et al.li2008t iKuntschn er et al.l 
l2oToh . In this work, we exploit the 4"2 diameter fibers and 
107 X 107" field of view of the Mitchell Spectrograph to study 
the spatial variation in age, metallicity, and abundance ratio 
gradients for eight massive early-type galaxies. Our increased 
leverage on stellar populations in the galaxy outskirts will allow 
us to put new constraints on the assembly of massive elliptical 
galaxies at late times. 

In §2 we describe the sample and in §3 we describe the in- 
strument and data reduction. The analysis is described in §4. 
Those most interested in results can focus on §5 and the sub- 
sequent discussion in §6. We summarize and conclude in §7. 
When needed, we us e the standard concordance cosmology of 
iDunklev et aT\i2QQ9ii . 

2. SAMPLE 

We start with a small pilot sample of eight galaxies as 
a proof of concept that the Mitchell Spectrograph is well- 
suited to this work (Table 1, Figure [T]l. The sample selec- 
tion is not ideal, and we do not make claims of its com- 
pleteness or uniformity, since the galaxies were selected with 
other science goals in mind. In short, we selec ted galaxies 
with red colors (u-r > 2.2; IStrateva et al.l 1200 ll) . stellar ve- 
locity dispersions that are larger than the instrumental resolu- 
tion of the Mitchell Spectrograph (ct*> 150 km s~'; as mea- 
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0.015 < z < 0.02 (85 Mpc, for a scale of 0.4 kpc per "). We 
examine all of the candidates and remove obvious edge-on disk 
galaxies, but have made no formal morphology cut. Thus there 
are two SOs in the final sample (CGCG137-019 and IC 1 153). 
We also make no selection on environment. However, in the 
sample there is a cluster galaxy (NGC 1270; iMiller & OwenI 
2001) , two brighte st group galaxies (NGC 677 and NGC 426; 
Berlind et al. 2006), two that belong to the same group (IC 11 52 
and IC 1153; White et al. 1999) and the rest are in lower density 
environments. Using our stellar velocity dispersion measure- 
ments, simple dynamical mass estimates for the galaxies range 
from 8 X 10"'-3 x 10" Mq, wit h only one belo w lO" Mq. 
Based on the mass function from Bell et alj (I2003h . they range 
from one to four times M* for ellipticals. 

3. OBSERVATIONS AND DATA REDUCTION 

The observations were obtained over two runs, one in Sept 

2010 (including the bulk of the galaxies) and the other in June 

2011 (IC 1152, IC 1153; Table 1). We used the George and 
Cynthia Mitchell Spectrograph (th e Mitchell Spectrograph, for- 
merly VIRUS-P; Imietal." 2008a) on the 2.7m Harlan J. Smith 
telescope at McDonald Observatory. The Mitchell Spectro- 
graph was built as a prototype for the VIRUS spectrograph that 
will soon be deployed on the Hobby-Eberly Telescope to per - 
form a dark energy experiment (HETDEX; iHill et al]|2008bl) . 
Each of the 246 fibers subtend s 4 "2 and are assembled in an ar- 
ray similar to Densepak (B arden et all 19981) with a 107" x 107" 
field of view and a one-third filling factor. The Mitchell Spec- 
trograph has performed a very successful search for Lyg emit- 
ters { Adams et al. 2011; Finkelstein et al. 2011; Bl anc et aTl 
l20Tlh and has become a highly productive tool to study spa- 
tially res olved kinematics and stellar populations in nearby 
galaxies iBlanc et al.ll2009l: lYoachimet allboiOl: iMurphv et atl 
l20lUlAdams et al.ll2012h . 

We used the low-resolution blue setting of the Mitchell Spec- 
trograph. Our wavelength range spans 3550-5850A with an av- 
erage spectral resolution of 5 A FWHM. This resolution deliv- 
ers a dispersion of ~ 1 . 1 A pixel"' and corresponds to (T» « 
150 km s~' at 4300 A, our bluest Lick index. Each galaxy 
was observed for a total of ~ 2 hours with one-third of the 
time spent at each of three dither positions to fill the field of 
view. Initial data reduction wa s accomplished using the cus- 
tom code Vaccine dAdams et al ]|2011:.Murphv et al..2011.) . We 
briefly review the steps of the pipeUne here, but refer the in- 
terested reader to the previous papers for more detailed dis- 
cussion. Initial overscan and bias subtraction are performed 
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Rest Wavelength (A) Rodius (R.) Rest Wavelength (A) Radius (RJ 

Fig. 1 . — Spectra (left, middle) and index equivalent widths (right) for our sample. Left, middle: Spectra are plotted in elliptical annuli of (from top to bottom) 
— 0. 5R(., 0.5 — li?c, 1 — 1.5i?c, \.5 — 2Re, and 2-2.5Re (see Table 1). Units are 10""erg s"' cm"^ A"', but the spectra have been offset for clarity. The sky feature 
in CGCG137-019 is highlighted with a dotted vertical line. Right: We show the Mgb index (green circles), the (Fe) index, which is the average of the Fe 5270 and 
Fe 5335 indices (red squares), and the H/3 index (blue triangles) in A. Error bars are derived via Monte Carlo simulations as described in §4.2.1. The (Fe) index is 
offset slightly in radius for clarity. 
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first on all science and calibration frames. All co-additions of 
data and calibration frames are performed with the biweight 
estimator (Beers et al. 1990). Twilight flats are used to con- 
struct a trace for each fiber, which takes into account curvature 
in the spatia l direction. W e employ a routine similar to that 
proposed by iKelsonI (l2003b to avoid interpolation during this 
step. Thanks to this special care, correlated errors are avoided 
and it is possible to track the S/N in each pixel through the 
remainder of the reductions. Knowing the S/N in each pixel en- 
ables deep er limits in detection experiments such as those per- 
formed bv lAdams et alj ( 1201 ll) . All subsequent operations are 
conducted in the new trace coordinate system within a cross- 
dispersion aperture of 5 pixels. 

To correct for curvature in the spectral direction, a wave- 
length solution is derived for each fiber based on arcs taken 
both at the start and end of the night. Gaussian fits to known arc 
lines are fit with a fourth-order polynomial to derive a complete 
wavelength solution for each fiber The typical rms residual 
variations about this best-fit fourth-order polynomial are 0.08 A 
for the Sept 2010 data and 0.04 A for our June 201 1 data. A he- 
liocentric correction is then calculated for each science frame. 

Next, a flat field is constructed from the twilight flats taken 
at both the start and end of the night. Variations in tempera- 
ture never exceeded 2 C for any of our observing nights and 
the stability of the flat fiel d has been shown to be < 0.1 pix- 
els under these conditions (lAdams et al.ll201 li) . As the twilight 
flats contain solar spectrum, we generate a mo del of this com - 
ponent by employing a bspline fitting routine ( lDierckxlll993h . 
A boxcar of 5 1 fibers is employed to model the solar spectrum 
that effectively removes all cosmic rays, continuum sources, 
and variations in the flat field, in order to isolate the solar spec- 
trum. The spatial and spectral curvature are leveraged here in 
order to supersample the solar spectra within the boxcar This 
supersampled bspline fit to the solar spectra within each fiber 
is then divided into the original flat. What remains are the flat 
field effects that we want to capture: variations in the individual 
pixel response, in the relative fiber-to-fiber variation, and in the 
cross-dispersion profile shape for every fiber. This flat field is 
then applied to all of the science frames. 

The next st ep is sky subtractio n. Unlike some instruments 
(e.g., Sauron; iBacon et"al] 1200 lb . the Mitchell Spectrograph 
does not have dedicated sky fibers. Instead, we observed off- 
galaxy sky frames with a sky-object-object-sky pattern, with 
five minute exposure times on sky and twenty minute object 
exposures. The sky nods are processed in the same manner as 
the science frames described above. In order to create a sky 
frame we give an equal weighting of two to each sky nod, then 
coadd them to achieve an equivalent exposure time as the sci- 
ence frames. The advantage of sky nods is the high S/N we 
achieve in our sky estimate, based on the large number of fibers 
in the IFU. The disadvantage is that we sample the sky at a dif- 
ferent time than the science frames and are thus subject to tem- 
poral variations in the night sky, particularly at twilight. How- 
ever, by varying the weights given to each sky nod we are able 
to explore possible systematics due to sky variability. For all 
lines measured, no EW value changes by more than 0.08 A due 
to temporal sky changes, even at the largest radii considered in 
this paper We give the details and return to possible systemat- 
ics related to sky subtraction in §4.2. 1 . Once the sky subtraction 
is complete, cosmic rays are identified and masked. 

We use softw are developed for the VENGA project 
jBlanc et alj|2009h for flux calibration and final processing. We 



observe flux calibration stars each night using a six-point dither 
pattern and derive a relative flux calibration in the standard way. 
Then we use tools developed by M. Song, et al. (in prepara- 
tion) to derive an absolute flux calibration relative to the SDSS 
imaging. M. Song uses synthetic photometry on each fiber and 
scales it to match the SDSS g-band image of each field, with 
a median final correction of ~ 20%. The correction exceeds 
50% only during a period of high cirrus in the second night 
of observing in Sept 2010, which affects NVSS J0320H-4136 
and NGC 426. Finally, all fibers are interpolated onto the same 
wavelength scale and combined. 

3.1. Radial bins 

We focus on spectra combined in elliptical annuli with a 
width of Q.SRe- Since the effective radii of these galaxies are 
~ twice the fiber diameter of 4", Q.5Re is roughly the scale 
of a single fiber In all cases we use the de Vaucouleurs ra- 
dius derived by the SDSS pipeline. Since we are averaging 
over such large physical areas on the sky, the exact measured 
Rg should not impact the conclusions. We did experiment with 
using wider radial bins at large radius to extend further from the 
galaxy center However, we did not boost the S/N appreciably 
in this way. Furthermore, we begin to be systematics limited at 

3Re (§4.2.1). The bins consist of 1-4 fibers at 0-0.57?,., in- 
creasing to 20-40 fibers at 2-2.57?^. The S/N per final coadded 
spectrum is shown in Table 2. 

Since all of our galaxies by selection have SDSS spectra, we 
can test the wavelength dependence of the flux calibration by 
comparing the shape of the spectrum in the central fiber of the 
MitcheU Spectrograph with the SDSS spectrum. We find ^ 5% 
agreement in nearly all cases, with no more than ^15% differ- 
ences at worst. 

4. ANALYSIS 

Our ultimate goal is to derive the stellar population prop- 
erties of the galaxies using the absorption line spectra. In 
principle, we can use fu ll spectral synthesis t e chniques (e.g., 
'Bruz ual & Charioi llool ICoelho et alJ l2007t IVazdekis et all 
2010), and exploit all the information available in the spectra. 
However, in practice, there are a number of hurdles, includ- 
ing imperfect flux calibration and systematic color effects, that 
make the model fitting sensitive in systematic ways to imperfec- 
tions in our data. It is known that the models are not always able 
to fit the absorption line equivalent widt hs (e.g ., Gallazzi e t al.l 
120051) . although see also iKoleva et al] ilOl lb . smce we are 
most interested in measuring changes in stellar population 
properties as a function of radius, we use Une index mea- 
surements (e.g.. Lick indices. iFaber et alii 19851 : IWorthev et al.l 
1994). As most prior work on this subject has used simi- 
lar methodology, we are in a good position to compare with 
the literature. Keep in mind that we are thus measuring the 
luminosity-weighted mean properties of the stellar population. 

Another major uncertainty in these stellar population synthe- 
sis models comes from our relative ignorance of stellar spec- 
tral energy distributions for stars with very different metal- 
licities and/or abundance ratio patterns than stars in our so- 
lar neighborhood. These modeling deficiencies impact our 
study directly, since elliptical galaxies tend to have higher 
metallicities and a-abundance ratios than local stars. Lately, 
the problem has garnered substantial atte ntion both from the 
point of view of full spectra synthesis (e.g., Coelho et al. '2007'; 
iConrov & van Dokkum .201 1; .Maraston & Stromback .201 1.) 
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Fig. 2. — a: Compaiison between the Lick indices measured from tlie SDSS spectrum and our central O.SRe bin for H/9 (blue triangles), the G-band at 4300A 
(purple squares), the Mgh triplet (green circles), the Fe 5270 line (red squares), and the Fe 5335 line (red hexagons). Fractionally, the two sets of indices agree 
within 3% of each other, while the mean deviation ((EWms — EWsdss)/EWsdss) = -0.02 ±0.08. The (Fe) index is systematically offset to lower values in our 
spectra by enough to introduce significant systematic offsets in the inferred metallicity. b: Fractional change in the EW of Ca H+K as a function of radius for all 
galaxies in the sample. Each color represents a different source. Note that the fluctuations between <Rf < 2 are within 10%, which is both the typical error in 
these measurements and the level at which Ca H+K is expected to vary due to real changes in stellar population. It is only in the outermost bin (2.5-3Re) that we see 
systematic effects begin to artificially lower the EW in some objects. The yellow points that rise at large radii coirespond to CGCG137— 019, which was observed at 
twilight and has impeii'ect sky subtraction. 



and for Lick index inversion meth ods jThomas et al. I I2003t 
ISchiavonll2007l: I Vazdekis et al.ll20T^ . We exploit these modern 
models, but note that our knowledge of the underlying stellai" 
evolution remains imperfect. 

4. 1 . Emission line contamination 

Something like 80% of elliptical galaxies contain low lev- 
els of io nized gas within an effective radius jSarzi et al. 
lYan &Blan ton 2011). For our purposes, this emission serves 
only as a contaminant, as it fills in the absorption features that 
we are trying to measure. By far the strongest emission fea- 
ture in our spectra is the [O II] AA3726, 3729 line, but there are 
no absorption features of interest that are confused by [O II]. It 
is contamination from H/? and [O III] A5007 that concerns us 
here. In order to correct for this low-level emissio n, we adapt 
the co de pPXF-nGANDALF developed by M. Sarzi ( Sar zi et al.1 
I2OO6I) and M. Cappellari (Cappellari & Emsellem 2004). pPXF 
performs a weighted fit to the galaxy continuum using spectral 
templates provided by the user, including a polynomial fit to the 
continuum and a Gaussian broadening to represent the intrin- 
sic dispersion of the galaxy. From these fits we derive a mea- 
surement of the stellar velocity dispersion in each radial bin. 
GANDALF iteratively measures the emission and continuum 
features simultaneously, to achieve an unb iased measurement 
of bot h components. As templates, we use iBruzual & Chariot! 
dlool single-age stellar population models with cr sa 70 km s 
resolution. In principle, we could use these fits to trace the stel- 
lar populations with radius, but for the reasons discussed above, 
we do not find this methodology robust. 

In general, the emission-line EW is small compared to that 
of the absorption lines. The maximum contamination comes in 
the central fiber of NGC 7509 (which has high-ionization lines 
indicative of an accreting black hole at the center). However, 
note that our spectral resolution of 150 km s"' is not necessar- 
ily high enough to resolve the emission Unes. 



[Mg/Fe] by 0.02 dex. Col. (9): [Fe/H] (dex). Col. (10|: [Mg/Fe] (dex). 

Table 2. Equivalent Widths and Inferred Stellar Population Properties 
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1.74±0.18 


427±0.24 


2,20±0.20 


2.16±0.38 


8.0 ± 1.70 


-0.43±0.13 


0.44±0.14 




2.5 


24 


1.33±0.22 


4.03±0.3I 


1.38±0.31 


l.S2±0.45 








NGC677 


0.5 


223 


1.57±0.02 


488±0.03 


2,56±0.03 


2.S9±0.04 


9,7 ± 0.20 


-0.20± 0.02 


0.38±0.02 




1.0 


139 


1.61 ±0.03 


457±0.05 


2,66±0.05 


2.62±0.08 


9,3 ± 0.30 


-0.1 9± 0.03 


0.30±0.02 




1.5 


76 


1.73±0.06 


426±0.09 


2.S5±0.08 


2.53±0.13 


7.9 ± 0.70 


-0.22± 0.06 


0.26±0.04 




2.0 


62 


1.62±0.08 


3.91±0.U 


2,44±0.11 


2.31 ±0.16 


9,7 ± 1.10 


-0.43± 0.07 


0.30±0.06 




2.5 


47 


1.53±0.10 


3.83±0.11 


2.43±0.14 


2J5±0.24 


11.9 ± 1.60 


-0.46± 0.10 


0.24±0.07 


NGC1270 


0.5 


189 


1.30±0.03 


5.38±0.03 


2,68±0.04 


2.S9±0.07 


15.0*± 0,50 


-0.27± 0.03 


0.46±0.02 




1.0 


145 


1.40±0.04 


5.31 ±0.06 


2,83±0.05 


2.76±0.07 


12.3*± 0.60 


-0.14± 0.03 


0.40±0.02 




1.5 


84 


1.46±0.05 


5.1 1 ±0.07 


2.42±0.07 


2.49±0.14 


ll.4*± 0.70 


-0.29± 0.05 


0.4S±0.04 




2.0 


57 


1.34±0.09 


4.92±0.09 


2,37±0.i0 


2.33±0.18 


15.2*± 1,00 


-0.45± 0.06 


0.50±0.0S 




2.5 


59 


1.33±0.10 


493±0.14 


2.67±0.I2 


2.49±0.20 


15.5*± 1.30 


-0.36± 0.08 


0.42±0.06 


NVSS 


0.5 


156 


1.50±0.03 


498±0.0S 


2.90±0.04 


2.76±0.06 


10.7 ± 0.50 


-0.1 2± 0.03 


0.32±0.02 




1.0 


39 


1.48±0.15 


460±0.16 


2,75±0.14 


2.31 ±0.24 


11.7 ± 1.90 


-0.42± 0.10 


0.44±0.10 




1.5 


47 


l.S0±0.13 


426±0.13 


2.25±0.13 


2.25±0.23 


6.6 ± 1.10 


-0.36± 0.14 


0.42±0.10 




2.0 


40 


1.45±0.12 


439±0.15 


2,06±0.16 


2.10±0.25 


13.5*± 1.50 


-0,58±0.12 


0.50±0.11 




2.5 


18 


1.77±0.35 


3.55±0.43 


1.92±0.41 


2.28±0.58 


7.3''± 2.30 


-0.41 ±0.24 


0.26±0.20 


IC1152 


0.5 


380 


1.75±0.01 


4.77±0.02 


2,95±0.01 


2.S6±0.03 


6.6 ±0.10 


0.03± 0.01 


0.26±0.01 




1.0 


129 


1.62±0.04 


4.4O±0.O6 


2,66±0.05 


2.61 ±0.08 


9,3 ± 0.40 


-0.20± 0.04 


0.26±0.02 




1.5 


60 


1.55±0.09 


417±0.13 


2.S0±0.I2 


2.41 ±0.1 8 


10.8 ± 1.20 


-0.38± 0.08 


0.30±0.06 




2.0 


49 


1.56±0.10 


436±0.15 


2,57±0.14 


2.S0±0.21 


10.6 ± 1.30 


-0.30± 0.09 


0.2S±0.0S 




2.5 


40 


1.43±0.12 


3.96±0.20 


2.4S±0.18 


2.45±0.30 


142*± 1.90 


-0.41 ±0.13 


0.22±0.09 


IC1153 


0.5 


536 


1.48±0.01 


418±0.02 


2,71 ±0.02 


2.72±0.03 


12.0 ± 0.20 


-0.23± 0.01 


0,16±0.00 




1.0 


128 


1.64±0.03 


3.99±0.04 


2.64±0.04 


2.65±0.05 


9.2 ± 0.40 


-0.20± 0.03 


0.16±0.02 




1.5 


71 


1.79±0.07 


3.92±0.n 


2.81 ±0.09 


2.68±0.16 


7.0 ± 1.00 


-0.13±O.O6 


0.14±0.04 




2.0 


61 


1.84±0.08 


3.78±0.1i 


2,93±0.11 


2.S4±0.16 


6,5 ± 0.90 


-0.02± 0.06 


0.06±0.04 




2.5 


44 


1.90±0.12 


3.57 ±0.16 


2.34±0.I4 


2.29±0.22 


5.9± 1.10 


-0.36± 0.10 


0.24±0.OS 


CGCGI37-019 


0.5 


151 


1.50±0.04 


4.26±0.05 


2,59±0.O+ 


2.47±0.06 


11.7 ±0.70 


-0.36± 0.03 


0.2S±0.02 




1.0 


94 


1.43±0.05 


415±0.07 


2.29±0.06 


2.28±0.10 


143*± 0.90 


-0.51 ±0.05 


0.34±0.04 




1.5 


52 


1.51±0.11 


3.71 ±0.13 


2.09±0.12 


2.13±0.19 


13.6*± 1.50 


-0.62± 0.08 


0.30±0.09 




2.0 


36 


1.45±0.15 


3.59±0.17 


2,08±O.I6 


2.ia±0.26 


15.3*± 2,10 


-0.60± 0.11 


0,22±0.09 




2.5 


24 


1.39±0.17 


3.89±0.26 


2.77±0.23 


2.47±0.38 


16.5*± 3.00 


-0.42± 0.16 


0.14±0.11 


NGC7509 


0.5 


187 


1.43±0.02 


4.84±0.O4 


2,44±0.03 


2.S0±0.05 


12,9*± 0,40 


-0.32± 0.02 


0,40±0.01 




1.0 


105 


1.45±0.05 


437±0.06 


2,48±0.05 


2.42±0.09 


13.0*± 0.60 


-0.40± 0.04 


0.34±0.O3 




1.5 


70 


1.39±0.07 


437±0.09 


2.37±0.08 


2.36±0.14 


148*± 1.00 


-0.47± 0.06 


0.36±0.06 




2.0 


53 


1.20±0.09 


3.98±0.13 


2,06±0.11 


2.32±0.18 










2.5 


43 


1.49±0.10 


3.57±0.]5 


1.83±0.15 


1.98±0.24 


13.9*± 1.60 


-0.76± 0.08 


0.42±0.08 


Note. — Col. (1 


Galaxy najne 


Col. (2): Radius (/;,.). Co 


. (31: Slgiial- 




ler pixel in [he c 


oiitiiiLiun at - 5 


OOOA for each 



radial bin. Col. (4): equivalent width (A), as measured by tick_ew. Tbe error bars are derived from Monte Carlo simulations and do not 
include systematic errors due to emission-line infill, etc. Col. (5): Mgb equivalent widtb (A). Col. (6): Fe 5250 equivalent width (A). Col. (7): 
(Fe) equivalent width (A). Col. (8): Age/(10' yr), as inferred by EZ^tges, based predominantly but not exclusively on the indices tabulated 
here. * indicates thai parameters were derived using H/3-I-0.2A. As explained in S4.3. we correct these ages by 33%, [Fe/H] by 0.09 dex, and 

We compare the fit to [O III] from the central pixel of the 
Mitchell spectra with the SDSS line fits at twice the resolution. 
The two fits agree within ^ 30% (range of 2-70% differences). 
As a result of the weakness of the features and the low spec- 
tral resolution, we do not have strong constraints on the flux 
or line shape of these emission lines. We incorporate emission 
line subtraction into our error bar estimates as described below. 
While the emission-line gas is not the focus of this paper, we 
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note that in a few cases the [O II] emission is observed to extend 
beyond IRe. An analysis of this emission will be the focus of a 
different work. 

Unfortunately, even very small corrections do lead to sub- 
stantive changes in the inferred galaxy ages - 0.1 A of H/3 
emission infill corresponds to an age difference of 1 Gyr 
([Schiavon 2007; Graves & Faber 2010). Thus, we caution that 
the absolute ages derived here are uncertain at this level. In 
some cases, in fact, the small infill of the H/? lines by emis- 
sion leads to H/3 EWs that are smaller than even the oldest 
and mo st metal-rich single stellar population (SSP) models. 
[Graves & Faber (2010) describes a detailed method to detect 
[O III] at the 0.2 A level. As described below, we do not have the 
S/N to perform their analysis, but the level of infill will become 
relevant again when we derive stellar populations in §4.3. 1. 

4.2. Equivalent Widths 

Measuring absorption-line indices and placing them on a 
common system is a delicate art that is sensitive not only to the 
spectral resol ution of the ins trument, but also flux calibratio n 
and S/N (e.g.. IWorfliev^tan [T994; Schiavon 2007: 1^120111) . 
We utiUze the flexible a nd robust IDE code lick_ew, written by 
G. Graves dGraves & Schiavon 2008). The code takes as input 
the stellar velocity dispersion at each radial bin and puts the 
measured EWs on the Lick system. In principle we measure 
EWs for all 26 Lick indices, but we focus our attention on H/3, 
Mgb, Fe 5270, and Fe 5335 (Table 2). We a lso m easure the 
Ca Hh-K index as defined by iBrodie & HanesI (Il986l) . Because 
of its high EW, any change in Ca Hh-K EW are indicative of 
systematic effects in the spectra (see §4.2.1). 

Since we compare with Lick indices defined on flux- 
calibrated stars dSchiavonl [2007 ). we only correct the indices 
to a standard spectral resolution but appl y no other zeropoin t 
offsets. This same approach is taken by [Graves et al.l (l2009h . 
By comparing the Mitchell indices from the central spectrum 
with the indices measured from the SDSS spectra, we confirm 
that we are on the same system (Figure |2]i. Note that the aper- 
tures are not perfectly matched (3" vs 4"2) but the difference 
between these two should be small since the observed gradi- 
ents are gentle. Spectra and Lick indices for all galaxies are 
presented in Figure [T] and Table 2. We find decent agreement 
between our indices and those from the SDSS spectra, except in 
the case of the Fe indices. Taking average differences AIndex = 
([IndexMitcheii-IndexsDSs]/IndexsDSs), and the standard devia- 
tion therein, we find AMgb= 0.03 ± 0.07, AH/3= -0.02 ± 0.09, 
AFe 5270= -0.07 ± 0.08, AFe 5335= -0.04 ± 0.08, and an 
overall offset of A Index= -0.02 ± 0.08 that also includes the 
G-band. 

If instead we look at absolute differences, we find that there 
is a small systematic difference in the Fe indices, in the sense 
that (Fcms) - (FesDSs) =-0.14A. While the systematic offset 
in (Fe) is small, it translates to large systematic errors of 0.08 
dex in [Fe/H] and 0.06 dex in [Mg/Fe]. We have explored var- 
ious causes for the systematic offset, including non-Gaussian 
line-broadening functions, variations in resolution with wave- 
length, and sky subtraction. None alone is sufficient to explain 
this systematic effect, although likely a combination of these, 
and possibly small-scale errors in flux-calibration, are to blame. 
As described above, our H/3 measurements are likely suffering 
from yery low levels of emission-line infill, which makes it dif- 
ficult to derive absolute ages. Thus, we do not focus on the 
absolute values of the derived parameters here, but rather on 



our main strength, which is gradients out to large radii. 
4.2.1. Uncertainties 

Potential contributors to our error budget include random 
noise, emission line removal, and sky subtraction. We can build 
the former two into Monte Carlo simulations of our measure- 
ment process. For each galaxy, in each radial bin, we start with 
the best-fit model from GANDALF and create 100 realizations, 
using the error spectrum generated by Vaccine. We rerun GAN- 
DALF and lick_ew on each artificial spectrum and take the er- 
ror bar as defined by the values encompassing 68% of the mock 
measurements. Note that the model contains both emission and 
absorption lines, and thus we should include the uncertainties 
due to emission line removal naturally in our error budget. 

Sky subtraction is a potentially significant source of uncer- 
tainty since we are working factors of a few below the level of 
the sky. In order to quantify how much sky variability effects 
our final science results we take a heuristic approach. We ex- 
plore two different scenarios. The first quantifies how changes 
in the night sky between our two sky nods impacts our final 
measured EW values. The second test quantifies the effect of 
an overall over- or under-subtraction of the night sky. In both 
cases we allow the weighting given to each sky nod to vary, 
then carry through all the resulting variations in the subtracted 
science frames. A comparison of the final measured EW values 
allows for a very direct measure of how well we are handling 
sky subtraction. We give the details of both scenarios here. 

To quantify how changes in the night sky between our two 
sky nods influence our final measured EW values we explore 
a range of weighting in the sky nods used for subtraction. For 
example, if the sky did not vary over the ^ 45 minutes between 
our two sky nods, an equal weighting of 2.0 given to each sky 
nod would be appropriate. This equates to an equal amount of 
total exposure time for each sky nod (e.g., 5 min x2.0 + 5 min 
x2.0 = 20 min). However, the sky may evolve on timescales 
shorter than this. To quantify the effect of an evolving sky on 
our final science results we ran several different sets of sky nod 
weightings through our reduction routines and compared the fi- 
nal EW values for all indices. We average deviations over radial 
bins from 1.5- 3Re- The sky weightings we explored varied by 
±20% from 2.0, yet always with a total weighting of 4.0 (e.g. 
1.8 for sky nod 1 and 2.2 for sky nod 2). We then make direct 
comparisons of the EW values (JEW = EWorig-EWnew) for all 
of the lines. As the observing conditions were predominantly 
stable for all our nights, we ran these tests on a single galaxy 
(NGC 677) and believe it to be representative. The largest SEW 
value measured was 0.080 ± 0.098 A (for H/3) when comparing 
a 2.2 - 1.8 to a 1.8 - 2.2 weighting. The SEW for each line for 
the case above were as follows: Mgb = 0.02 ±0.02 A: H/3 = 
0.080 ± 0. 1 A: Fe 5270 = -0.007 ± 0.02 A. From an analysis of 
th e variability of the sk y spectra over several nights (Figure 15 
m iMurphv etal.ll20T T) the case of 20% variability is extreme. 
The one exception occurred with IC 1152, which saw a rising 
moon for some of the exposures. In this case we explored a 
wide range of sky weights and found a weighting of 2.4 - 1.6 to 
be optimal. 

The second scenario we explore is aimed at understanding 
what systematic effect over- or under-subtraction of the sky has 
on our EW values. To test this we conducted a similar set of 
tests, yet allowed the final weighting of 4.0 to vary. We ran 
tests exploring both a 5% and 10% over and under subtraction. 



Elliptical Outskirts 



1 




-0.6 -0.4 -0.2 0.2 0.4 0.6 



[Fe/H],3 [a/Fe],3 

Fig. 3. — We compare m etallicities ([Fe/H]; a) and a abundances ([a/Fe] as measured by the Mgb index; b) derived using the prescriptions of lGraves & SchiavonI 
120081) and IThomas et al J j2p03) respectively. We show observations for all galaxies at radii out to 2.5Re. Each galaxy is represented by a different color. The 
agreement between the two models is decent, with ([Fe/H]Qs— [Fc/HItmb) = —0.06 ±0.1, and ([a/Fe]Gs-[a/Fe]TMB) = —0.03 ±0.08. To reduce crowding, we 
include error bars for only two galaxies in the sample. 



relative to equal exposure time. We then made the same com- 
parison in (5EW as described above. In the case of a ±5% sys- 
tematic error in subtraction, we find 5Mgb of -0.09 ± 0.04 A, 
of 0.05 ±0.09 A, and 5<¥s> of -0.05 ±0.02 A. The worst 
deviations in individual bins are at the 0.15 A level in M.gb and 
H/3 for the 5% oversubtraction case. It is interesting to note that 
there are not strong systematic effects. Instead we see the in- 
dices bounce around at the 0. 1 - 0. 15 A level for this level of sky 
subtraction error. When we get to 10% over-subtraction, the 
errors are 5Ugb of -0.2 ± 0.08 A, SnjS of 0.01 ± 0.04 A, and 
(5<Fe> of -0.04 ±0.08 A. The worst deviations in individual 
bins are at the 0.2 A level in Mg/? and H/3 for 10% oversubtrac- 
tion. Larger fractional errors in sky level would lead to obvious 
residuals in our outer fibers that we do not see. As all of the 
5EW values calculated from both tests described here are within 
our typical uncertainties, and the scenarios we tested were ex- 
treme cases, we conclude that our results are robust against sky 
variations on the scale of ^ 45 minutes seen in our data set. 

As an additional sanity check of our sky subtraction, we 
calculate the EW of the Ca Hh-K AA3934,3968 Hnes. These 
features have very high EW, but also are virtually insensi- 
tive to changes in ste llar populations (at the ^ 10% level; 
iBrodie & HanesI (Il986l) ). Thus, we expect the line depths to 
be constant with radius. Because these features are quite blue, 
they provide a rather stringent test of our fidelity in the sense 
that the blue spectral shape is most sensitive to errors in sky 
subtraction. We find that these lines do not vary by more than 
10% out to 2.5Re for most systems (FigurelU?). As a result, we 
view our results out to 2.5Re as reliable. While we show the 
point at 2.5-37?^ in the figures, we will not use that point in our 
fitting. 



4.3. Stellar population modeling 

We now convert the observed EWs into ages, metallicities, 
and abundance ratios using SSP models. Since all indices are 
blends of multiple elements and all depend on age, meta llic- 
ity, a nd abundance ratio to some degree (e.g., Worthev e t al.l 
1 19941) . modeling is required to invert the observed EWs and 
infer stellar population properties. We compare two differ- 
ent mode ling techniques. W e first use the methodology out- 
lined in iThomas et al.l (|2003, 2005). Since the EWs of Mgb 
and (Fe) change with both [Fe/H] and [Mg/Fe], these authors 
construct linear combinations of the two indices. The index 
[MgFe']= vTMg/7[0.72Fe5270 + 0.28Fe5335]) is independent 
of [Mg/Fe] and tracks [Fe/H], while the index M.gbl(Pt) de- 
pends only on [Mg/Fe]. The pair of derived indices is then in- 
verted to infer metal content and ab undance ratios. 

We also use the code EZ_Ages bv lGraves & SchiavonI (l2008l 
see also Schiavon 2007). Here, the age, metallicity, and al- 
pha abundance ratios are fit iteratively using the full suite of 
Lick indices. EZ_Ages solves for the best-fit parameters by 
taking pairs of measured quantities (e.g., H/? EW and (Fe)) 
and then locating the measurements in a grid of model val- 
ues spanning the full range of age and (in this case) [Fe/H] 
abundance of the models. The model has a hierarchy of 
measurement pairs that it considers, first pinning down age 
and [Fe/H], then looking at [Mg/Fe] and so on. The code 
then iterates to improve the best fit values. Many more el- 
emental abundances can be fitted by EZ_Ages, thus enabling 
study of the independent variability of [N/Fe], [C/Fe], etc. 
There is evidence that individual a elemental r atios vary inde- 
pendently in individual Milky Way stars (e.g., 'Fulbriaht et al.l 
2007) and possibly i n galaxies as well (e.g.. Kelson et al.j 
l200a ISchiavonI |2007|) . We do not have adequate S/N in 
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Fig. 4. — Radial profiles of Age (10'" yr; blue triangles), [Fe/H] (red squares), and [Mg/Fe] (green circles) for each galaxy. Radius is measured in units of Re. 
Open symbols are the adjusted values using H/3+0.2A, and the two gaps are cases that still fell off of the grid. As above, error bars are derived from Monte Carlo 
simulations. 
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the blue in dices to derive other elemental abundances (e.g., 
lYanI I2Q III) . po we i ust i assume tlfat [Mg/Feil tracks [q /Fel . 




log (Model^B+oa/Model) 



HB+o.a/ 

Fig. 5. — The difference between model pai'ameters from EZ_Ages from the 
observed data and when the H/3 index is boosted by 0.2 A. When H/3 is artifi- 
cially increased, the ages (blue solid histogram) are lower by 31 it 5%, [Fe/H] 
(red dashed histogram) increases by 0.09 ± 0.01 dex and [Mg/Fe] (green long- 
dashed histogram) increases by 0.02 ± 0.01 . When the H/3 EW is too low and 
falls off of the grid, these corrections are used to put the model parameters 
based on H/3-I-0.2A on the same scale. 



For the remainder of the paper we will use [a/Fe] to refer to the 
a-abundance ratios collectively, bearing in mind that we have 
directly measured [Mg/Fe]. 

The Thomas et al. (TMB) approach and the Graves & Schi- 
avon (GS) model have somewhat different philosophies, but are 
inherently similar Both are based on the inversion of single- 
burst model grids. Both adjust their primary models for variable 
g-abundanc e ratios at a range of metallicities. T MB use solar 
isochrones ( Cassisi et al.ll997tlBono etal ■1199 7') but then mod- 
ify the indices using the response functions of Tripicco & Belli 
(Il995h . GS use solar isochro nes from Girardi et al. (2000) and 



Qf-enhanced isochron es from [Sa lasnich e t al., (,2000,) . with the 
response functions of iKorn eta l. (2005). TMB invert a small 
number of high S/N indices, while GS rely on all measured in- 
dices in an iterative fashion. In principle we can test some of 
the systematics of the modeling by comparing our results from 

the two different approaches^ 

Finally, note that Graves & SchiavonI (|2008|) parametrize 
their models in terms of [Fe/H] rather than total metallicity 
[Z/H], so that th ey report direct observables. As described in 
ISchiavonI (l2007l and references therein), we are not able to 
measure [O/H] directly, and thus cannot truly constrain [Z/H]. 
To c ompare the tw o models, we will use the standard conver- 
sion dTantalo et al.iri998UThomas et al.ll2003h : 



[Fe/H] = [Z/H]-0.94[a/Fe] 



(1) 



In Figure [3] we compare [Fe/H] from Thomas et al. 
([Fe/H]TMB) with that from Graves & Schiavon ([Fe/H]Gs) and 



likewise for [Mg/Fe], using observations in all radial bins for 
each galaxy. The agreement is reasonable. Overall, we see a 
scatter of ( [Fe/H] gs- [Fe/H] tmb) = -0.06 ±0.1 for [Fe/H] and 
([Q!/Fe]Gs-[a/Fe]TMB) = -0.02 ±0.08 for [a/Fe]. For the rest 
of the paper we will focus on the EZ_Ages results from Graves 
& Schiavon, but trust that our results can be directly compared 
with many in the literature. Also, we have rerun the EZ_Ages 
modeling with the solar isochrones, and the derived metallici- 
ties and a-abundance ratios agree within the measurement er- 
rors. The radial dependence of the derived quantities for each 
galaxy is shown in Figure]?] 

4.3. 1 . Low Hj3 Equivalent Widths 

As mentioned above, there is very likely real but undetectable 
levels of H/3 emission that slightl y lowers the obse rved H/3 
EWfl At the levels measured by iGraves & Fabed (1201 Ol) in 
composite SPSS spectra (~ 0.2 A), the age errors are ~ 2 Gyr 
in general (ISchiavonI ]2007]) . In some cases, we cannot derive 
reasonable model parameters because the measured H/3 EW is 
too low to fall onto the SSP grids. Since we do not have the S/N 
needed to correct our spectra on a case by case basis, we use the 
following procedure to derive model parameters at the radial 
bins where the H/3 index fall off the bottom of the grid. Note 
that these corrections are not strictly correct, since the level of 
emission must vary with radial distance. However, they are the 
best that we can do at present. 

We recalculate the age, [Fe/H], and [Mg/Fe] for all galaxies 
with the H/3 EW increased by 0.2 A. We then derive an average 
difference in each measured property between the two sets of 
models, as shown in Figure ]5] The differences in [Mg/Fe] and 
[Fe/H] are very small and (crucially) show no trend with radius, 
S/N, or H/3 index. The run with increased H/3 EW returns ages 
that are 31 ± 5% dex lower, [Fe/H] values that are 0.09 ± 0.01 
dex higher and [a/Fe] values that are 0.02 ± 0.01 dex higher on 
average than the unadjusted data. 

We correct the model parameters derived from the H/S+Q.IA 
run to align with the fiducial models using the corrections listed 
above. At each radial bin where we could not derive model pa- 
rameters using the fiducial H/3 EWs, we instead utilize the cor- 
rected model parameters from the H/3-1-0.2 A run. In the follow- 
ing, all such points are indicated with open, rather than filled, 
symbols. With this correction, we derive SSP properties for 
all but two radials bins over all of the galaxies. Again, apply- 
ing the same H/3 correction for all points that fall of the grid 
is not strictly correct, since there is likely radial dependence in 
the amount of infill. Therefore, we again emphasize that our 
main strength is in measuring the radial trends rather than the 
absolute values of the stellar population parameters. Note that 
Kelson et al. (2006) take a similar, although perhaps more nu- 
anced, approach by shifting all of their model grids to match 
the Lick indices of their oldest galaxies. 

5. THE AGES AND METAL CONTENTS OF STELLAR HALOS 

The most striking trend in Figure JT] is the clear and steady 
decline in the Mgb index out to large radii. While the gradients 
vary from object to object, the qualitative behavior is the same 
for all systems. In contrast, both the (Fe) index and the H/3 in- 
dex are generally consistent with remaining flat over the entire 



^We have also investigated whether a non-Gaussian line-broadening function or changing spectral resolution could lead to H/3 infill. We broaden the Bruzual & 
Chariot models with the appropriate line-broadening function, and find that the measured indices only change at the ~ 0.05 A level. As described above, in principle 
sky subtraction could cause errors at the 0.1 A level, but it is hard to understand how those eiTors would be so systematic. We conclude that low-level emission is the 
most likely culprit. 
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radial range. In Table 3 we show the gradients in Hf3, (Fe), and 
Mgb measured asdX = S log X/S log R/R^ for each index "X". 
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Fig. 6. — Correlation between the Mgb index and central stellar velocity 
dispersion. For direct comparison with the literature, we use the stellar veloc- 
ity dispersion measured within an effective radius, although the differences in 
dispersion are small as a function of radius (see §5.1 below). We plot measure- 
ments for each galaxy at three radii: the central O.SRe (small red filled circles), 
between 1.5 and 2 Re (medium yellow circles), and between 2 and 2.5 Re (large 
blue circles). The SO galaxies CGCG137-()19 and IC 1153 are indicated with 
a double cir cle on the central (re d) point. For reference, we show the measure- 
ments from lGraves et al.l 120091) as small dots. Since these are SDSS galaxies 
at z 0.1, the hght comes from ~ Re. We also show the average relation from 
iTrager et alj I200ffl as dashed lines, which show the trend line well within the 
effective radii of average elliptical galaxies. In terms of Mgh EW, the halos of 
these elliptical galaxies have a similar chemical makeup as galaxies of much 
lower mass. 

We now ask whe ther the well-known Mgb-a^. relation (e.g., 
iBender et al.ll 19931) . is preserved at large radii. In Figure|6]we 
plot the Mgb index at O.SRe, l.SRe, and l.SRe as a function 
of galaxy stellar velocity dispersion measured within the effec- 
tive radius. This figure visually displays two interesting trends. 
First of all, the Mgb EWs beyond IR,, in these massive ellipti- 
cal galaxies fall significantly below the central Mg^j-cr* relation. 
Matching the Mgb EWs at IR,, with the centers of smaller ellip- 
tical galaxies suggests that the halo stars were formed in smaller 
systems. In §5.2 we find that if these stars were accreted from 
smaller elliptical galaxies, they would come in ^ 10 : 1 merg- 
ers. Of course, the more detailed abundance patterns of the halo 
stars will give us more clues as to the possible origins of these 
halo stars. 

Secondly, we do see hints of a Mgb-a^, correlation even at 
large radii, but with much more scatter There is also an intrigu- 
ing hint that the slope of the Mgb-cr^, relation changes. Unfortu- 
nately, the correlation is driven to a large degree by the galaxy 
NGC 1270, which has the largest cr* value in our sample. NGC 
1270 also happens to be our only cluster galaxy. Since the dis- 
tribution of merger mass ratios is nearly independent of mass 
jFakhouri et al. 2010), we might expect the mass of the typical 
accreted system to rise with mass, thus preserving an Mgb-a^, 
relation at large radius. However, we will need a larger sample 
at the highest velocity dispersions to say for certain whether a 



Mgb-a^, trend continues in the galaxy outskirts. 

The next obvious question is whether the Mgb EW drops pri- 
marily because of changes in [Fe/H] or [a/Fe]. As outlined in 
the introduction, all evidence suggests that metallicity decline 
is the primary cause for the decline in Mgb EW. Since very little 
data exists at such large radii, however, it is worth investigating 
the [Fe/H] and [a/Fe] measurements directly. Again we empha- 
size that the absolute values of the derived [Fe/H] and [a/Fe] 
are unreliable, but that the gradients should be robust. Below 
we present gradients in the derived metallicities and abundance 
ratios in order to determine what drives the striking decline in 
Mgb EW. 

5.1. Gradients 



Table 3. Lick Index, Age, Metallicity, and Abundance Ratio Gradients 



Galaxy 


log a. 


6 HP 




5 Mgb 


S (Fe) 






|Fe/H| 


S [Mg/Fe] 


U) 


(2) 


(3) 




("4) 


(5) 










NGC426 


2.44 


-().10±0.21 


-0 


15±0.02 


-0.09±0.06 


0.25±0.17 


-0 


30±0.I6 


-0.30±0.13 


NGC677 


2.37 


O.OSiO.IO 


-0 


ISiO.OI 


-0.08±0.03 


-0.1 3 ±0.09 


-0 


24±0.07 


-0.26±0.06 


NGC 1270 


2.49 


0.1I±0.I2 


-0 


lOiO.OI 


-0.06±0.03 


-0.20±0.13 


-0 


I6±0.08 


-0.12±0.06 


NVSS 


2.37 


0.!0±0.2I 


-0 


16±0.02 


-0.23±0.05 


-0.22±0.17 


-0 


75 ±0.15 


0.26±0.12 


ICII52 


2.36 


-().!5±0.13 


-0 


15±0.0I 


-0.17±0.03 


0.64 ±0.10 


-0 


67±0.08 


-0.09±0.06 


ICII53 


2.40 


0.20±0.09 


-0 


12±0.0I 


-0.07±0.02 


-0.76±0.09 





16±0.06 


-0.19±0.12 


CGCG 137-019 


2.23 


-().08±0.18 


-0 


13±0.02 


-0.14±0.04 


0.40±0.23 


-0 


47±0.11 


-0.06±0.09 


NGC7509 


2.28 


-().07±0. 1 1 


-0 


18±0.0I 


-0.11 ±0.03 


O.I7±O.I5 


-0 


49±0.08 


-0.14±0.07 



Note. — Col. fl): Galaxy name. Col. (2): Log stellar velocity dispersion within the effective radius (km .s ). Col. (3): 
Logarithmic gradient in W index 5H3/-51og (ff/ff^). Col. (4): LogariUimic gradient in Mgb index (5Mgb/(Slog (fl/ff,.). Co!. 
(5): Logarithmic abundance gradient 5(Fe)/51og (J?/flc). Col. (6): Logarithmic age gradient 51og (age/Gyr)/51og tR/R^). Col. 
(7): Logarithmic metallicity gradient 5[Fe/H]/51og (ff/Rr). Col. (8): Logarithmic abundance gradient 5[Mg/Fe]/5iog (fl/flf ). 



As described above, we measure the gradients 6Age, 
(5[Fe/H], and (5[Mg/Fe] as SX = 6 log X/S log R/R^. We use 
adjusted values of age, [Fe/H], and [Mg/Fe] for objects that fell 
off of the grid due to low Hj3 EW. We do a very simple least- 
squares fit to the measured quantities as a function of logarith- 
mic radii in units of 7?^ (Table 3). Only [Fe/H] shows signifi- 
cant evidence for a significant radial gradient. On average, the 
[a/Fe] ratio declines very gently with radius, but the trend is 
not significant in individual cases. The age gradients take both 
positive and negative values, but are rarely significant. Given 
the uncertainties with H/3 described above, we will focus ex- 
clusively on metallicity and abundance ratio gradients here. 

We should note that unlike most observations in the litera- 
ture, our fits are weighted towards the outer parts of the galax- 
ies, and we have less spatial resolution in the galaxy centers. 
We are thus less susceptible to stellar-population variations in 
the central regions caused by late-time accretion and/or star 
formation. Indeed, iBaes et al.l (l2007l) find a clear break in the 
slopes of metallicity gradients at smal l radii, with the gradi ents 
getting shallower at larger radii, as do lCoccato et al.l (1201 Oh for 
the brightest cluster galaxy in Coma, NGC 4889. Our mild 
metallicity gradients are similar to those seen at large radius by 
these authors, as well as by G. Graves & J. Murphy in prepara- 
tion in M 87. 

In principle, correlations between the stellar population gra- 
dients and other properties of the galaxy provide additional 
clues as to the origin of the gradients. We examine the relation 
between g^, and gradients in F igurelT] Our results are consistent 
with the lSpolaor et al.l (1201 Ol) result that the gradients in ellipti- 
cal galaxies show a larger scatter at higher cr, . In our small sam- 
ple, we do not find a correlation between the decline in metallic- 
ity and the isophote shape, with the metallicity gradients being 
substantially shallower than the decline in isophote level with 
radius. We do not find a correlation between gradients in stel- 
lar velocity dispersion and gradients in indices or gradients in 
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metallicity or abundance patterns, but the sample is yet small. 
Eventually it will be interesting to look for correlations between 
the local escape velocity and the metallicity gradient as may be 
seen if the meta llicity is set locally by the ability of gas to es- 
cape the galaxy (IScott et al.ll2009l;IWeiimans et al.ll2009l) . 



5.2. Where do the Halo Stars Come From? 

Figure |6] presents the intriguing possibility that we may un- 
cover the mass scale of accreted satellites by matching the 
metallicities and abundance ratios of the stellar halos with the 
centers of smaller elliptical galaxies. In this section we take 
the measured gradients in Mgb index, [Fe/H], and [Mg/Fe] and 
attempt to constrain the typical mass of an accreted satellite. 
As we will show, the abundance patterns in the halo stars do 
not match the central regions of any local elliptical galaxies, 
thereby complicating our efforts to find the progenitors of the 
halo stars. Nevertheless, we can derive an approximate mass 
ratio from our measured gradients, accepting that present-day 
galaxies do not form perfect analogs of the accreted satellites. 

Using our observed ct* within R,,, we assign a central value 
of abundance ratio and metallicity using the Mgb-(7t, [F e/H]- 
cr*, and [Mg/Fe]-(T* relations from iGraves et alj (l2007l) . We 
then use our observed gradients to link the stars at large radii 
with the (7* of its most likely progenitor. For instance, we as- 
sign an [Fe/H] to ea ch galaxy using the [Fe/H]-(T* relation from 
[Graves et al.l (l2007h . Then we use our measured gradients to 
calculate [Fe/H] at 2-2.5Re- That [Fe/H] value at large ra- 
dius is matched to the cr* of an accreted galaxy with the same 
metallicity using the [Fe/H] -cr* relation. The derived value of 
(7* is translated into a stellar mass M* using the projections 
of the Fundamental Plane presented by Desroches et al. (2007). 
The typical mass of accreted galaxies based on each property is 
shown in Figure [8] 

Figure [8] quantifies the typical mass of accreted satellites. 
The figure strengthens our conclusion from Figure |6] Based 
on the Mgb EWs alone, we find that the stars at > 2Re were 
accreted from galaxies ^10 times less massive than our target 
galaxies (blue points). However, the figure makes very clear 
that the metallicities and abundance ratios of our stellar halos 
cannot be matched with present day elliptical galaxies of any 
mass. Looked at another way, it says that the [Fe/H] gradients 
we see in individual halos are steeper than the [Fe/H] -cr* rela- 
tion for the population overall, while the [Mg/Fe] gradients are 
shallower (Figure |7]). The Mgb-a^, relation in the galaxy out- 
skirts taken alone suggests that stellar halos are built by ^ 10 : 1 
mergers, but the halo stars have lower metallicities and higher 
a-abundance ratios than present-day low-mass ellipticals. 

Of course, accreted satellite systems need not be small ellip- 
tical galaxies, and may well have started their lives as gas-rich 
disk galaxies. Thus, it is useful to consider spiral galaxies as 
well. Since disks also have low surface brightness, there are few 
studies of ab undance patterns in disks based on stellar absorp - 
tion features jGanda et al.ll2007t lYoachim & Dalc antonir2008l) . 
Furthermore, the interpretation using SSP models is severely 
complicated by the clear o ngoing star formation in these sys- 
tems ( Mac Arthur et alj|2009i) . There is general agreement that 
massive red bulges show simila r patterns and scaling rela- 



tions as elliptical g alaxies (e.g., Moorthv & Holtzman 2006t 
iRobaina et"al]|201lh . At lower mass, Ganda et al. (2007) find 
that even the bulge regions of later-type spirals have younger 
ages, lower metallicities and solar abundance ratios compared 
to elliptical galaxies. 



Due to their protracted star formation histories, present-day 
late-type spirals do not share the high abundance ratios of the 
stellar halos studied here. Instead, it may be more productive 
to consider individual components of present-day galaxies. For 
instance, thick disks are clearly older than thin disks, and de- 
pending on their formation channel may share characteristics 
of these stellar halos. However, aside from our own galaxy, 
it remains very difficult to obtain robust stellar abundance pat- 
terns in thick disks (.Yoachim & Dalcanton 2008.) . In Figure 
|9]we show the abundance ratios at 2-2.5Re in our sample as 
derived using the central relations from Graves et al. and our 
gradients. We compare with other stellar populations including 
elliptical and spiral galaxies, and subcomponents of our own 
galaxy. This figure shows that some of our galaxy halos have 
similar abundance patterns as the Milky Way thick disk while 
others are more consistent with low-mass elliptical galaxies. 
We now address plausible scenarios for how the stellar halos 
were assembled. 

6. DISCUSSION 
6.1. Theoretical Expectation 

We have established clear gradients in Mgb EW with radius 
in all of the galaxies in our study. Based on Lick index inver- 
sion methods we have argued that in general these gradients are 
dominated by metallicity with a weak contribution from abun- 
dance ratio gradients as well. While the Mgb EWs of the galaxy 
halos match those of galaxies an order of magnitude less mas- 
sive, the stars appear to have lower metallicity and higher a- 
abundance ratios than do low-mass ellipticals. We now review 
the various physical processes that we believe can impact the 
observed chemistry, in order to determine which scenarios are 
favored by our observations. 

We start with "monolithic" collapse, by which some large 
fraction of the galax y is built in a single dissipational burst of 
star formation (e.g., Egg en et"aL|[T962l) . While we believe that 
we live in a hierarchical Universe in which large galaxies are 
built up through the merging of smaller parts, the central stel- 
lar populations of massive ellipticals clearly i mply that their 
stars were formed rapidly at redshifts z>2 (e.g..l Thomas et all 
l2005i and references therein). A rapid dissipational phase at 
high redshift is likely an important part of elliptic al galaxy 
formation, followed by late-time dry merging (e.g. , iTal et al.l 
'2009t lvan Dokkum et allbOlOtlNewman et al]|201 ih . There is 
strong theoretica l support for such a "two-phase" picture (e.g., 
iNaabet all [20091: lOser et alJl20Toh and high-redshift progeni- 
tors, in whi ch star formation h as ceased at early times, are ob- 
served (e.g.. lKriek et"Zll2009l) . 

A large body of work, starting with Larson] d 19741) . has con- 
sidered the chemical evolution of monolithic collapse models, 
including a heuristic star-formation law, chemical enrichment, 
and (typically) galaxy-scale mass loss driven by supernovae. 
Since these galaxies have deep potential wells, the gas in the 
center cannot easily be ejected, but instead is enriched by pre- 
ceding generations of star formation and grows metal-rich. In 
contrast, in the outer parts of the galaxy, winds can be effec- 
tive at ejecting metals. Steep gradients in metallicity and a 
abundance ensue (e.g., JOirlberg 1984; Arimo to & Yoshiill 19871: 
iKawata & Gibsonll20 03T Kobavashi 2004). 

In a modern cosmological context, massive halos 
that host the progenitors of massive elliptical galax- 
ies are constantly bomba rded with smaller halos (e.g., 
iBovlan-Kolchin et all l2009h . Furthermore, we see galaxies 
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Fig. 7. — (a): The Mgb index as a function of the stellar velocity dispersion measured at each radial annulus rather than within Re. Each galaxy is shown as 
a different colored circle. The congelation between ct, and Mg b in individual galax ies as a function of radius is considerably steeper than that seen across galaxy 
centers. As in Figure[6]above, we show the measurements from 'Graves et alj (2009D as small dots. Since these are SDSS galaxies at z 0.1, the light comes from 
~ Re. We also show the average relation from Trager et al. (2000), which shows the trend line well within the effective radii of average elliptical galaxies, (b): 
Gradients in metallicity 5[Z/H]/<5 log (R/Re) as compared with the stellar velocity dispersion of the galaxy. Our galaxies are shown as large red circles. We compare 
with the compilation of Spolaor et al. ( 2010), and find reasonable agreement in the velocity dispersion range covered by our data. We measure cr, within Re, while 
Spolaor et al. use 1 /SRe, but the gradients in ct» over these radii are generally small (e.g., Jorgensen et al. 1995; Cappellari et al. 2006). Also note that Spolaor et al. 
tabulate gradients out to the effective radii in most cases, while we measure them between O.SRe and 2.5Re. (c): Same as (b) for gradients in (5[o/Fe]M log (R/Re) 
(big green circles). Again note that [Mg/Fe] is assumed to trace [a/Fe]. 



merging dToomre & Toomr^l 19721: ISchweize3ll982h . Mergers 
will scramble the orbits of the constituent parts to some degree 
and wash out metallicity gradients (e.g., I White! [1980 ). The 
degree of mixing depends on a variety of factors. In a violent 
relaxation scenario (without dissipation) existing stars may not 
migrate much, th us preserving the original chemical patterns 
jvan Alb ada 1982). In contrast, gas-rich major merging should 
efficiently supply gas to the center of the remnant, where it will 
form metal-enriched stars and steepen metall icity gradients 
(e.g.. lMihos & Hernquisll996l:ICox et al.ll2006 >). Modern sim- 
ulations that account for cosmological merging and chemical 
evolution conclude that merger remnants will have shallower 
metallicity gradients on average, wit h a much larger scatter 
than the monolithic collapse case (e.g.. lKobava"shill2004 . 

6.2. Constraints From Data 

Let us briefly review what the observations tell us (see also 
the introduction for more complete references): 

1 . There is a strong correlation between M.gb EW and ct* 
measured w i thin the effective radius of the galaxy (e.g . , 
iFabed 119731: iDressler et al.lll987t iBender et al.lll993h . 
The majority of this trend is attributed to metallicity, 
but ther e is also a trend between g-abundan ce ratio and 
cr* (e.g.. lWordiev et al.1 19921: [Graves et al.f20 09). Much 
like the mass-metallicity trend observed in star forming 
galaxies, these trends must arise at some level because 
of the relative ease of ejecting metals f rom the shal- 
low potentials of low-mass galaxies (e.g.. lDekel & SiU3 
119861: iTremonti et alll2004l) . 

2. In ^ L* elliptical galaxies, the metallicity decreases 
outwards gently, falling by ^0.1 -0.3 dex per decade 
in radius. The gradients are too shallow in general 



to agree with pure monolithic collapse scenarios. No 
clear trends are seen in a-abundance ratio gradients, 
with increasing, decreasing and flat trends observed 
(e.g., iKuntschner et al.l 1201 Oh . again in disagreement 
with monolithic collapse. In more massive elliptical 
galaxies, metallicity gradients still dominate, but there 
is a wider di spersion in the gradient slopes at a given 
(T* dCaroUo et al. 1993: Spolaor et al.„2010i) . consistent 
with what we observe (Figure |7). 

3 . We add an additional robust spectroscopic point beyond 
2Re. The Mg/j EW in the stellar halos match the values 
seen in local elliptical galaxies that are ten times less 
massive than our sample galaxies. However, the centers 
of present-day elliptical galaxies in this mass range have 
more metals and lower values of [a/Fe] than do the stel- 
lar halos. Our gradient observations are strongly in con- 
trast with predictions from monolithic collapse scenar- 
ios, in which g-a bundances would decrease outwards 
(lKobavashill2004l) . The question is whether we can ex- 
plain the observed stellar properties if the halos were 
built via minor merging at late times. 

At first glance, our observed a-abundance ratios are difficult 
to understand in any scenario. We rule out a pure monolithic- 
collapse scenario because of the lack of gradient in abundance 
ratios. Major mergers are not strictly excluded, but seem un- 
able to produce such consistent decreasing metallicity gradi- 
ents without some tuning. If, in contrast, the outskirts are built 
up via minor merging at z < 1 (e.g., Naab et al. 2009), then 
we would expect the stars to have the same metallicities and 
abundance patterns as small elUptical galaxies today. Figure |8] 
demonstrates that [a/Fe] is too high at a given [Fe/H] to derive 
from present-day low-mass elliptical galaxies. 
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Fig. 8. — The characteristic mass of an accreted satelhte at 2.5Re as inferred 
from the Mgh EW (blue circles), [Fe/H] (red squares), and [Mg/Fe] (green tri- 
angles). As in Figure|6] SOs are indicated with double circles. To derive the ac- 
creted mass, we assume prev iously derived Mgh-Ut., [Fe/H]-(T, , and [Mg/Fe]- 
<jf relations IGraves et al.l2007,) to assign central values of these quantities to 
our galaxies. Using our observed gradients, we match the abundance patterns 
beyond 2Re to a present-day elliptical galaxy with the same Mgfc EW, [Fe/H], 
or [Mg/Fe] . We translate th e er« value to a st ellar mass using the Fundamental 
Plane relations from Desro ches et alj j2007l) . The total galaxy mass derived 
using the Fundamental Plane is shown with the solid black line. We fit av- 
erage relations between cr* and Ma(.(.,etcd based on each measurement with a 
fixed slope to guide the eye only (Mgfc dashed blue line; [Fe/H], dot-dashed 
red line; [Mg/Fe], long-dashed green line). Our indirect technique allows us to 
circumvent uncertainties in the absolute values of [Fe/H] and [Mg/Fe] . If the 
stellar halos were constructed from analogs of present-day ellipticals, then the 
accreted mass derived from each indicator would agree. Instead we find that 
present-day ellipticals cannot simultaneously match the observed low values 
of [Fe/H] and high values of [Mg/Fe] . 

It is useful to draw an analogy with studies of the Milky Way 
halo. Early suggestions that the halo may be built by the accre- 
tion of satellites (Searle & Zinn 1978) were called into ques- 
tion by the observation that the abundance patterns of stars in 
the Milky Way halo do no t match those of the existing satel- 
lites (e.g. JVenn et"al]|2004 and references therein). However, 
satellites that are accreted early by the Milky Way will have 
a truncated star-formation history. Thus, they will have high 
[a/Fe] ratios compared to systems that continue to accrete gas 
and form stars u nt il the present day (e. g . , [Robertson et al.l2005l: 
iFont et al.ll2006l) . iTissera et all (1201 ll) track the chemical evo- 
lution of eight Milky-Way analogs with hydrodynamical sim- 
ulations. Most of the mass in each simulated halo is accreted 
(much of it after a redshift of one) rather than formed in situ. 
The accreted stars are metal poor and a-enhanced, since the 
stars were predominantly formed at early times in small halos 
that have truncated star formation histories. We observe the 
same trend in the halos of the massive early-type galaxies ex- 
amined here. 

Our galaxies are considerably more massive than the Milky 
Way, and for the most part they do not have disks at the present 
time. There must be many differences in their merger histories 
from that of the Milky Way, and the absolute metallicities of the 
Milky Way halo stars are much lower (Fig. |9]l. Nevertheless, 



we suggest that the halos in both cases are built by the accretion 
of smaller galaxies before these small systems have a chance to 
self-enrich. Thus, we support a scenario in which massive ellip- 
tical galaxi es were built up v ia gas-free minor merging at late 
times (e.g.. iNaab et al.ll2009l) . The abundance patterns of the 
massive halos differ from those seen in ^ L* ellipticals today 
because the former had their star formation history truncated 
when they were accreted. In contrast, we know that ^ L* el- 
lipticals had some ongoing star formation at late times (e .g., 
iBabul & Ree"i[T992llT^omas et al.ll2005HKoleva et al.ll201 Ih . 



The simulations presented in lOser et all (|201dl see also 
C. Lackner et al. in preparation) provide strong support for halo 
build-up via the late accretion of small satellites. They consider 
galaxies with M* ranging from 5x 10"' - 4x 10^' Mq , and find 
that half of the stellar mass was accreted after a redshift of z « 1 . 
They are able to reproduce the observed size evolution in ellip- 
tical galaxies (e.g., van der Wei et al. 2008). Furthermore, they 
reproduce the observed increase in stellar mass on the red se - 
quence over the last eight billion years (e.g.. lFaber et al.ll2007l) . 
Of more direct importance to our story, while accreted late, the 
majority of the accreted stars were formed at z > 3. These stars 
were by necessity formed rapidly out of low-metallicity gas. Of 
course, our sample is still small, but eventually we hope to have 
a large enough sample to look for differences in elliptical halo 
abundance patterns as a function of mass. 

At a redshift of z w 1, the progenitors of local L* ellipticals 
could well have consisted of old, metal-poor, and a-en hanced 
stars. In the schematic picture of iThomas et al. I (l2005l) . ~ L* 
ellipticals are forming the majority of their stars over a few Gyr 
around z < 1 . If roughly half of the stars are formed at z < 1 , 
with a timescale of ^ 3 Gyr, while the original population has 
[Q!/Fe]=0.2-0.3, the final mass-weighted abundance would be 
[a/Fe]sa0.15, as observed forL* ellipticals today. Thus, those 
galaxies that are not cannibalized by larger galaxies can eas- 
ily self-enrich to form the observed populations today. We 
make a testable prediction of strong evolution in the metallicity 
and [a/Fe] ratios of ^ L* elliptical galaxies from z ~ 1 to the 
present. 

In our analysis, we have assumed a constant initial mass 
function. If instead massive elliptic al gal axies have either a 
top-heavy (van Dokkum 2008; Dav^ l2008l) or a bottom-heavy 
dvan Dokkum & Conrovn201 1 ) initial mass function compared 
to lower-mass systems, that may change the interpretation of 
the gradients. We have also ignored the possible variation 
between different a elements, particularly nitrogen CSchiavonI 
l2007h . 

7. SUMMARY 

We have used the Mitchell Spectrograph to gather high S/N 
spectra of eight massive early-type galaxies out to 2.5Re, in or- 
der to study the chemistry of their stellar halos. Looking first 
at the trends in Lick indices with radius, we find that the EW 
of Mgb drops, such that the well-known Mg/j-cr* relation is not 
preserved at large radii. Instead the Mgb EWs at large radii are 
similar to those found in the centers of galaxies that are an order 
of magnitude less massive. 

We show that the well-known metallicity gradients seen 
within Re continue to the furthest radii probed here. In con- 
trast, [a/Fe] does not drop substantially in any object, and cer- 
tainly never approaches the solar value. Thus, the stars in the 
outer regions of these elliptical galaxies are metal-poor and a- 
enhanced, much like the stars in the Milky Way halo. We sug- 
gest that the outer parts of these galaxies are built up via minor 
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merging with a ratio of ~ 10 : 1, but that the accreted galaxies 
did not have sufficient time to lower their a-abundance ratios to 
those seen in ^ L* elliptical galaxies today. 
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Fig. 9. — We i nfer the [Fe/H] and [Mg/Fe] values at 2—2.5Re using the cen- 
tral relations of IGraves et al.ll20Q7h and our measured gradients (large blue 
circles). We compare the a bundance ratios an d metallicites in our stellar halos 
with Milky Way stars fromlVe nn et alj 120041) , including thin disk (small black 
circles), thick disk (small grey open squares), and ha lo (small grey stars) stars. 
For comparison we also show the track of the 'Graves et al. ( 2003) composite 
elliptical galaxies from the SDSS (filled red squares) and the central regions of 
late-type spiral bulges from Ganda et al. 1,2007, filled blue triangles). Again, 
taken as a group, our stellar halos are not well-matched by the integrated prop- 
erties of galaxy centers today. However, we do see some overlap with Milky 
Way thick disk stars and low-mass elliptical galaxies. 



This paper is only a proof of concept; the Mitchell Spectro- 
graph is ideally suited to study the faint outer parts of galax- 
ies, and there is a considerable amount of follow-up work to be 
done. First of all, we would like to investigate the kinematics in 
the outer parts of these galaxies to determine whether there are 
correlations between angular momentum content and metallic- 
ity. We are working on gathering a larger sample, with a full 
sampling of velocity dispersion, size, and environment, to see 
whether the radial gradients in (e.g.,) metallicity, correlate with 
the size of the galaxy at fixed cr*, or the large-scale environmen- 
tal density. It seems clear that galaxy e volution is accelerated 
in ric h environments (e.g.. lThomas et a l. 2005; Papov ich et al.l 
1201 lb . leavi ng subtle impr ints in the central stellar populations 
of galaxies dZhu et alj|201 0). Whether that will leave clear sig- 
natures in the gradients remains unknown. Even in our own 
small sample there are real differences from galaxy to galaxy, 
and it will be very interesting to see whether the large-scale en- 
vironment is the cause. 

The referee gave us an extremely prompt and thorough re- 
port that considerably improved this manuscript. We thank G. 
Blanc and M. Song for crucial assistance with data reduction. 
We thank G. Graves, J. E. Gunn, L. C. Ho, J. P. Ostriker, and 
B. E. Robertson for many stimulating discussions about both 
the measurements and the science. J.M.C. is supported by an 
NSF Astronomy and Astrophysics Postdoctoral Fellowship un- 
der award AST- 11 02525. 
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